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I. Governing Equations 



Elastic nod 



i 



T = -I 



lTii(xi - t) 

^ 




Area A 



'C=TiTii(xi + -2i,t) 




Xi — 2" ^1 ^1 + ~2" 



■> 6(xi, t) 



Thin elastic rod with axis In the Xx-dinection and uniform cross secdon of area A 
(a] ttie rod; (b) force and tractions applied to an element of length Ax^ centered 
at Xi with displacement 5(X]^, t). 



5^5 

AAx.p^ 



Til Xi + — ^ -TiJ X. 



AXi 



A + F A AXi 



d^5 _ a Til 

et^ ~ 5xi 



+ F ( F is body force density) 



55 



T^i = E Hooke's Law (stress-strain relation) 

Young's modulus (modulus of elasticity) 



5^5 ^ 5^5 ^ 
= E^ + F 

dt^ dxj 
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Modulus of Elasticity E and Density p for Representative Materials 



Material 



E-unlts of lO'-^N/m^ p-units of lO^kg/m-' Vp-units+ m/sec 



Aluminum (pure and alloy) 


0,68-0,79 


2.66-2.89 


5100 


Brass (60-70% Cu, 40-30% Zn) 


1.0-1.1 


8.36-8.51 


3500 


Copper 


1.17-1.24 


S.9E-8.9S 


3700 


Iron, cast (Z. 7-3.5% C) 


0.89-1.45 


6.96-7.35 


4000 


Steel (carbon and iow alloy) 


1.93-2.20 


7.73-7.87 


5100 


Stainless steel (18% Cr, 8% Ni) 


1.93-2.06 


7.65-7.93 


5100 


Titanium (pure and alloy) 


1.06-1.14 


4.52 


4900 


Glass 


0.49-0.79 


2.38-3.88 


4500 


Methyl methacrylate 


0.024-0.O34 


1.16 


1600 


Polyethylene 


1.38-3.8 X 10"^ 


0.915 


530 


Rubber 


0.79-4.1 X 10-= 


0.99-1.245 


46 



Sec S.H Crandall, and N.C. Dahl, An Introduction to W?e Mechanics ofSo!idi=, McGraw Hill, New York, 1959, for a 
list of references for these constants and a list of these constants in English units. 
-|- Computed from average values of E and p. 



II. Example Case Study - Change in Rod Length 



^1 



/ 



Cross-sectional 
area A=lcm^ 



l=0.1m 



ATiiCI) 



f=100 lb=445 N 



Metal rod fixed at Xi=0 and subject to the force f at Xi=l 
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steady state: — = 0, F = 

5t 



g^5 



=> 5 = axi + b 



Boundary Condition at x^^ = 1 : AT]^i(xi=l) = AE 



as 



= AEa =f 



Xi.| 



5(xi = 0) = = b 



5(xi) 



fx 
AE 



5(xi=l) = 



AE 



E = 0.7x10^^ N/m^ (Aluminum) 



A = 0.0001 m^, f = 445 N, I = 0.1m 



5(xi=0 



445 (.1) 



6.36xl0^m = 6.36 ^m 



10-4 (.7x10^^) 



III. Elastic Waves : F = 



A. Wave Equation 



5^5 E 5^5 



et^ P ex 



2 ' 'P 



V = E/ 
P V/P 



5 = 5^ (x-Vpt) + 5_ (x + Vpt) 
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Superposition : 5 = 5+ (x - Vpt) = 5+ (a) , a = x - Vpt 



95_95 da _ ^ 55 
9t 9a 9t 9a 



9^5 
9t2 



9a 



8 8] da 



P 9aj 9t 



2 9^5 



P 9a2 



95_95 9a_95 
9x 9a dx da 



9^5 
9x2 



9 (95 



9 f9519a 9^5 



dx\daj da\da) dx da^ 



9^5 



2 9^5 2 9^5 2 9^5 

= V = V = V 

9t2 P9a2 ^9x2 P9a2 



Q.E.D 



For 5 = 5_ (x + Vpt) let Vp -Vp 



B. Velocity and Stress Variables 



95 ^ ^95 
V = — , T = E — 
9t 9x 



9^5 _ 9v _ 9T 
9t2 ~ 9t ~ 9x 



9T 
9t 



9^5 



9v 



9x9t 9x 



_9_ 
9t 



_9_ 
9x 



' 9t 



9T 
9x 



9T _ g 9v 
9t ~ 9x 



9^7 
9x9t 



92v 



= E 



9^/ 
9x2 



92v 



2 S^V 2 c/ 
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d_ 
dx 



5t 



1 dT 
p 5x 



1 d^T 



1 d^T 



E St' 



at 



5v _ 1 5T 

~ E "sT 



5^ 



/' ^ v' 

p ex' ' 



E/P 



V = (x - Vpt) + v_ (x + Vpt) 
T = T^x-Vpt) + T (x + Vpt) 



C. Method of Characteristics 



Lagrange 
l^ultipliers 



ev r,ev ^ dT dT ^ 
p — + — - — = 

at ax at ax 



av ^av aT „ aT 
— - E — + — + — = 
at ax at at 



^av_^av ^ 
-^^ at p>.i ax 



dv 
dt 



^ar Xi aT 
p^i at 8x 



dv(x,t) = — dt + — dx 

^ at ax 



ciT(x,t) = — dt + — dx 
^ at ax 



dx 
dt 



-X2 E 



E 
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—7 + — = ^ dv ± -= dT = 

jit pX-^ ^pE 



dv + -= dT = on 



dv — ;= dT = on 



v + ^^T = C 



on 



dx 
dt 

dx 
dt 

dx 

d^^-^P 



on 



dx 

dt =""P 



V = 



C + C 



T C - 



D. "Cook-Book" Method 



J / ..\ 5v . Sv Sv ST 

dv(x,t) = — dx + — dt , p — = — 
^ ' 8x dt dt 8x 



.-rf ^^ ST. 5T.^ 5T_5v 
dT(x,t) = — dx + — dt , — =E — 
^ ^ 8x dt dt dx 



P 
-E 



-1 

1 



dt dx 
dt dx 





'5v/et " 




"0 " 




5v/5x 









5T/5t 




dv 




dT/dx 




dT 
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- 1' 









-E 


1 







aet 














dv 


dx 












dT 





dt 


dx 


dt 




p 








- 1 









-E 


1 







det 














dt 


dx 


















dt 


dx 



set each determinant to 
zero for general solution 



1^' Characteristic equation 



det 



p 0-1 

-E 1 

dt dx 

dt dx 



-E 
dx 




1 









dt dx 



+ dt 




-E 1 





dt dx 



= pdx 



"-E 1" 


- dt 


~-E 1' 


dx 




dt 



-(dxf p + E(dtf 



'dx Y _ E dx 
.dtJ " p"^ dt 



±/e7^ = ±v 
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2"'^ Characteristic equation 



det 



"0 





-1" 




-E 


1 







dv dx 










dT 


dt 


dx 




^ Tdv " 
dT dt 


+ 


"dv 
dT 


dx" 




-E 1 
dv dx 
dT dt 



= Edvdt -dx dT = 



dx E 
Edv - dT — = ^ Edv+ j- dT = 
dt "V P 



dv + = 



V + 



constant 



on 



dx 
dt 



+ v„ 



V + -;= = C_ on 
VpE 



dx 
dt 



IV. Case Study : Region of Initial Unifornn Stress 



t = 0: T(x,t = 0) = 



Tm X < a 



X > a 
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AW2 
< 



^ 





Thin rod sub)ect to an Initial uniform static stress 1,^ over the section -a<x<a. At t=0 the externa I 
forces ATn,/2 are removed. 
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VIII. Boundary Conditions 




13 



AT(l,t) 

i 



AT(l,t)+Bv(l,t)=0 

C-= C+(B-A>fpE ) 
B+A\fpE 



Matched end: no 

reflection 

B=AsrpE 



BvCI,t) 



Simple boundary conditions on the end of a thin rod: (a) free end; (b) fixed end; 
(c) end attaclied to a damper producing a totai force Bv. 
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t 



-I — L 



Boundary where 
T{l, t)=0 




t 



Free end of rod 



Boundary where 
T(-l, t)=0 



1^ 



(a) Thin rod of length 21 centered at x=0; (b) an x-t plot showing the characteristics 
relevant to the effect of the ends on the dynamics. 
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Propagation of an Initial pulse of stress on a rod terminated at x = ± I with free ends. 



IX. Sinusoidal Steady State 



A. Velocity Driven Thin Rod 



v(x=0, t) = VoSlncct >■ 







Rxed end 



v(x,t) = Re v(x)e 



B.C: v(x = 0) = - jvo 
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v(x = l) = 



5^v 2 S^v 



= V 



-co^vCx) = 



dx 



2 "P 



d^V i2"/' \ r>i2 2/2 PtC^ 

— ^ + k^v(x) =0, =«) /Vo - ^ 
dx2 V ^ ' / P 



v(x) = Asinkx + Bcoskx 

v(x = 0) = B = -jVg 

v(x = l) = 0= AsinkI - JVq coski => A=jVoCOtkl 



v(x) = i Vn I ^^5-!^ sinkx - coskx 
^ ^ ' sinki 



= ■^^^ (coskI sinkx -coskx sinkl) 
sinkI ^ ^ 

jVoSink(x-l) 
sinkl 



5T i \ y c 1 cosk(x-l) 
— = E — => ia)T(x) = j Vo E k V; — ^ 



. ^ k cosk(x-l) 

T(X) = Vn E i ^ 

^ ^ ° CO sinkl 



= VqEVp/E 



cosk(x - 1) 
sinkl 



= VqVpE 



cosk(x - 1) 
sinkl 



v(x,t) = Re [v(x)e 



. . I sinkfx - I) 

I = - Vo 5^ ^ sinrat 

sinkl 
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T(x,t) = Re[f(x) 



Jcot I ^ V\/Pi (-Qs|</x-l)c0S(0t 

sinki ^ ' 



Resonance: sinkI = => kl = nn => k = 



(0 _ nn 



X. Electromechanical Coupling 




s 2 ds 



1 v^SqA, 



fx has steady part & transient part 



1 Vq' Sq a, 

2 (g.5(-l,t))^ 



5 = 5o (x) + 5 (x t) 



g 



steady part I ^§0 = 



1, 25(-l,t) 

g 
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= ^ Sq = ax + b 



5o (x = 0) = b = 



= ax 



' dx 



= + 
x=-i 2 



1 SqA.Vo^ 



2al 



EAr a 



EA, 



So A, yi I 



1 So A, Vo^ 

+ — - — s— ^ => a 



+ -£0 



EAr - 



So A, Vo^l 



Transient Part: 8' (x,t) = Rep(x)ej'"'] 



5' 



5' (x) = A;^ si n kx + A2 cos l<x, l< = co/Vp 



5(x = 0) = = A2 



6' (x) = A^ sinl<x 



B.C. at X 



EA, 



dx 



^^^|^5'(-i,t) 



x=-l 
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So A, Vo" I 



< I ^ k = jk| 



tanki =tan(jkil) = jtanh(kil) =jkil 
(Oi = k,v„ 



A I- « 3 A 

EA,g^ ^ 



Unstable toot 




eiO<ait) = e-«Ht (Unbounded if e)j<0) 



— tanhlql 
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